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We employ Bloch wavefunctions similar to those used in calculating electronic and optical 
band gaps in solid state physics to derive dispersion relations for water wave propagation 
in the presence of an infinite array of periodically arranged surface scatterers. For one 
dimensional periodicity (stripes) , band gaps for wavevectors in the direction of periodicity 
are found corresponding to certain frequency ranges which support only nonpropagating 
standing waves as a consequence of multiple Bragg scattering. The dependence of these 
band gaps as a function of scatterer density, strength, and water depth is analyzed. 
We find that in contrast to band gap behavior in electronic, photonic, and acoustic 
systems, these gaps can increase with excitation frequency to. Thus, higher order Bragg 
scattering can play a dominant role in suppressing wave propagation. Furthermore, in one 
dimension, an additional constraint (in addition to single scatterer energy and momentum 
conservation) on the calculation of transmission T and reflection R coefficients of a finite 
number of scatterers is calculated from the exact dispersion relation of the infinite periodic 
system. This relationship may be useful in measurements where only phase or amplitude 
is easily measured. In simple two dimensional periodic geometries, we find no complete 
band gaps, implying that there are always certain directions which support propagating 
waves; we offer semi-quantitative reasons for why this is so. The role of evanescent modes 
is discussed, and finally, equations for water wave band structure in the presence of a 
uniform fiow are derived. 



1. Introduction 



Wave diffraction is ubiquitous in Nature and important for countless technological 
applications. One application is sea surface wave interactions with natural and man- 
made structures. Coastlines as well as off-shore structures are susceptible to wave 
forces, and much research has focused on diminishing these interactions with the aim 
of designing more effective off-shore rigs, breakwaters, and docks (see for example 
Mei, Hara fc Naciri (1988)], |Mei (1983)| and [Porter fc Evans (1995)|). Wave interactions 



with ice in polar waters is another important aspect of surface wave scattering. Different 
ice floe structures and distributions modify the sea surface wave spectra and must be 
considered when analyzing remote and direct wave sensing data and models of wind wave 
generation and oceanic transport ( Wadhams et al. (1986) ). 

The theories behind wave diffraction have a rich history of development in many areas 
of physical science. In most scattering problems, for example those involving acoustic, 
electromagnetic, or electronic wave functions, the governing equation is a second order 
linear wave equation of the Helmholtz type. Usually, two boundary conditions are satis- 
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fied at the interface of scattering bodies, along with a radiation condition far away from 
the body. These prescriptions uniquely solve most scattering problems. 

However, although surface perturbations on an ideal fluid propagate as waves, and 
general principles of scattering theory from other contexts apply, the equation govern- 
ing the fluid velocity in the bulk liquid has no spatial variation. Thus, scattering can 
occur only through interactions with localized geometrical variation of the boundaries, 
or localized variations in the conditions imposed at the boundaries. An example of 
the former is t he scattering of water waves from variations in bottom depth studied 



by |Mei (1985)|, [Belzons, Guazzelh fc Parodi (1988)| , [Davies, Guazzefli fc Belzons (1989) 



and Mitra fc Greenberg (1984) . This problem has many i mportant applications. 



including underwat er breakwaters (Mei, Hara fc Naciri (1988) ), long distance ocean 



wave propag ation (Elter fc Molyneaux (1972) ) 



near beaches ( |Mei (1985)| , [Naciri fc Mei (1 



wave propagation over rippled b eds 
Davies, Guazzelli fc Belzons (1989)| and 



Kirby (1986) ), and third sound propagation of superfluid helium films as discussed by 



Kleinert (1990) . Characterizing surface waves can also be a means of remote sensing 



of bottom topographies. For abrupt jumps in bottom depth geometry, eigenfunction 
expansions of velocity and pressure are usually matched at these regions. The infinite 
number of poles in surface water wave dispersion is an indication that near scatterers, an 
infinite number of evanescent modes are required to match boundary conditions (see Eqn 
(Al)). Periodic bottom depths, or periodically arranged impenetrable obstacles piercing 
the sea surface, have also been studied theoretically with multiple scattering analyses 
by |Mei (1983)| , [Porter fc Evans (1995] , [Kagemoto fc Yue (198*6]. Integral equation tech- 
niques have been applied to study periodic breakwaters (Fernyhough fc Evans (1995)) 
and bottom depths (Black, Mei & Bray (1971)). In these studies the possibil- 
ity of resonant Bragg scattering from finite structural arrays was recognized by 

[Mitra fc Greenberg (1984) ; however, 



Mei (1983), Davies, Guazzelh fc Belzons (1989) 



only first order Bragg scattering has been discussed. The problem of random bot- 
toms can also be studied, especially using shallow water theory, where interesting 
loc alization behavior has been theoreti cally predicted and experimentally obse rved 
by Devillard, Dunlop fc Souillard (1988) , and Belzons, Guazzelli fc Parodi (1988) and 



Elter fc Molyneaux (1972) , respectively. 



In this paper we consider scattering that occurs when the parameters in the bound- 
ary conditions are spatially varying along the interface. Physical examples can be 
found on widely varying scales where surface tension or surface bending rigidity 
modulate periodically. For example, Lucassen-Reynders fc Lucassen (1969) show that 
surfactant deposition at the air-water interface drastically attenuates wave propaga- 
tion, and multiple surface wave scattering from surfactant concentration variations is 



thought to be one contributing mechanism as suggested by Chou and Nelson (1994) 
Moreover, monolayers of surfactants at air-water interfaces can phase separate into 
periodic domains as a consequence of line tension and microscopic dipolar interac- 
tions (Andelman, Brochard fc Joanny (1987) and Knober & Desai (1992)). Similarly, 
ice floes, which vary the surface mass density as well as contribute bending moments 
to the air-sea interface, also scatter water waves and attenuate wave energy in Marginal 
Ice Zones (MIZ) ( Fox fc Squire (1990) and Squire et al. (1995) ). Here, regular patterns 
of surface ice thickness can form due to periodic convection of freezing surface waters, 
or current instabilities of a freezing mushy ice zone. Also, stress from incident water 
waves on an ice sheet can break strips of ice off the leading edge and form very long, 
regular stripes of ice separated by open sea ( Squire et al. (1995)1 ). An approximate two- 
dimensional periodicity also exists in regions of pancake ice, where circular floes occupy 
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(a) 

Figure 1. Schematic of periodic scatterers at the Uquid interface. The depth is a constant h, 
and the bending rigidity and/or surface tension is D2 and/or (T2 inside the hatched regions and 
Di and/or ai outside. Depicted here are scatterers which have sharp discontinuities in D and 
a along the surface, in (a) one- dimensional stripes and (b) square array of disks of radius Rq. 



a high filling fraction of the surface. All of the physical systems mentioned are schemat- 
ically depicted in Figure |l| in the limit of perfect periodicity. 

The con nection to the classic p robl em of single o bstacle scattering studied for ex- 
ample by Greene fc Heins (1953) and Ursell (1957) is worth mentioning. The limit 
of thin surface obstacles, the prototypical problem of a semi-infinite or finite dock, 
has been studied extensively u sing various sophisticated applied mathematical meth- 
ods such as integral equations (Friedrichs fc Lcwy (1948) ) and Wiener-Hopf techniques 



( [Greene fc Heins (1953) , Heins (1956) ). In these problems, the fiuid normal velocity at 
a simple surface region (half-space or strip) is constrained to vanish. Surface bend- 
ing moments however, have been shown to be important for the study of water wave 
propagation in flexible docks and ice by Liu fc MoUo-Christensen (1988)| . Reflection and 
transmission coefficients of water wave scattering from a semi-infinite flexible dock have 
been calculated using variational methods and numerical solution of integral equations 
by Squire et al. (1995) and Meylan fc Squire (1994) ; however, only the amplitudes have 
been explicitly calculated. When considering multiple scattering, however, the phase is 
a critical variable to be considered as it affects whether waves interfere constructively 
or destructively. Wiener-Hopf techniques applied directly to the differential equations 
(Noble (1958), Heins (1956)) of motion are inapplicable when a semi-infinite dock has 
surface tension or bending rigidity since the proper solution requires a higher order dis- 
continuity at the boundary between free water and the dock edge to conserve vertical 
stresses or bending moments, and ultimately energy (this failure is apparent when one 
considers the discontinuity in surface tension (t{x) giving a d{x) in Eqn. (^.8[)). The 
more rigorous integral formulation of the Wiener-Hopf method by Krein (1962) may 



prove useful in solving the semi-infinite, flexible dock problem. An integral equation 
solution to the reflection and transmission coef ficients of scattering from a sem i-inflnite 
surface tension discontinuity has been found by |Gou, Messiter fc Schulze (1993)| , but the 
energy conservation was applied a posteori and phase shifts were not found. Meylan 
fc Squire 1994 have calculated the diffracted wave field from a circular sheet of flexible 
ice by matching many evanescent modes in a Green's function expansion. The circular 



Stone (1995) 



surface tension contrast problem has also been treated by Chou, Lucas 
using dual integral equations. 

Despite the technical difficulties in calculating scattering from single scatterers, we will 
show that an exact solution to an infinite periodic array of surface scatterers (such as ice 
floes) can be simply calculated. The results represent the effects of multiple scattering 
and allow easy physical interpretation. When used with certain robust wave properties 
such as energy and momentum conservation, and radiation conditions, the solutions can 
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yield useful relationships in the single surface scatterer problem. Numerous other physical 
results, interpretations, and implications for applications can also be gleaned. 

In the following Section we formulate the model describing surface wave propagation 
in the presence of an infinite periodic array of surface deformable scatterers (Fig. We 
adapt standard techniques used in solid state physics to describe periodic electronic wave 
functions to study surface water waves. Other applications of periodic scattering arise in 
the study of acoustic and photonic band gaps in periodic inclusions and colloidal crystals, 
respectively. In our problem, the surface boundary conditions supply the periodicity and 
are cast into matrix equations with the dispersion relation of an infinitely periodic system 
determined by the spectrum of eigenvalues. 

In the Results and Discussion, Section 3, we consider both one and two-dimensional 
arrays of surface scatterers; wave propagation of course, can occur in any direction in the 
(x, ?/)-plane (see Fig |l|). The physical features of the dispersion relations are displayed 
and qualitatively explained. We find b and gaps, (also often denoted as stop gaps) analo- 



gous to those found in othe r systems ( Ashcroft fc Mermin (1976) , McCall et al. (1991) 



Bigalas fc Economu (1992) ), where as a consequence of the coupling between the peri- 
odic nature of the scatterers and the plane waves, the dispersion curves w(g) have jumps 
between which the excitation frequencies do not correspond to extended plane waves, 
but rather localized states corresponding to complex wave vectors. The location of these 
jumps in a; as a function of wavevector q define Bragg points (ID) or lines (2D) in 
Fourier space where resonant multiple scattering occurs (Ashcroft & Mermin (1976) and 
Joanopoulos, Meade & Winn (1995)), and large distance propagating modes cancel. For 
one-dimensional periodicities, an approximate form is given for the splitting of doubly 
degenerate eigenvalues into a gap. We emphasize the physical phenomena of band gaps 
and their dependences on parameters such as scattering strength and water depth, we 
will not provide an exhaustive array of numerical calculations. We also briefly discuss 
the differences between water wave scattering and wave scattering derived from a "local" 
Helmholtz equation. In particular the effects of the decaying "evanescent" components 
near the scatterers are discussed. 

In the Summary and Conclusions, we briefly mention the general features of the ba- 
sic results and their implications for multiple scattering in natural and manmade ap- 
plications, such as wave propagation in ice fields and arrays of surface plates. In the 
Appendices, we use a remarkably simple but general analysis of one dimensional sym- 
metric scatterers, to derive in certain limits, the band gaps as a function of transmis- 
sion/reflection coefficients and their phases. Therefore, upon comparison with the exact 
results of Section 3, we obtain an additional constraint on transmission and reflection co- 
efficients (including their phases) of a finite number n one-dimensional symmetric surface 
scatterers of arbitrary strength. We also outline the equations required to treat period- 
ically surface scattered water wave propagation in the presence of underlying uniform 
flows. 



2. Model of Periodic Surface Waves 

In this section, we formulate the the dynamics of an interface overlying an ideal fluid 
of depth h (Fig. Balance of surface stresses leads to dynamical equations which 
govern surface wave propagation. We restrict our analysis to surface scatterers. These 
structures are assumed to have small thicknesses d (d <^ X, h, where A is any wavelength) 
and lie entirely at the water interface z ~ and thus do not change the domain over 
which Laplace's equation for the velocity potential holds. Examples are thin plates with 
bending stiffness or domains of surface active materials which locally decrease surface 
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tension. The scatterers are periodic and their in-plane interfacial positions are assumed 
fixed. 

2.1. Surface Water Waves 

The motion of the interface is coupled to the bulk flow which is described by the linearized 
Navier-stokes equation for incompressible ideal flows, 

9tv + (v • V)v = --Vp + pg. z<ri(f,t) (2.1) 
P 

where p is the bulk fluid density, pg = —pgz is the external body force density due 
to gravity, p is the pressure, and r]{f,t) is the surface deformation as a function of 
surface coordinate f= {x,y). We neglect all dynamical effects of the upper fluid, air. 
By assuming irrotational flows, the velocity can be expressed in terms of a potential, 
V = Wip. Hence, 

V-v = VV(^^,i) = (Vi+9^^)¥'(r,2,i) = 0. z < Tj{r) (2.2) 

where = + dy is the 2D Laplacian in the coordinates of Fig. |l|. The linearized 
kinematic conditions at the free interface and the bottom z ~ —h are 



and 



\im dz(p{f,z,t) ^ dtri (2.3) 



lini ^ d^(p{f, z, t) = 0. (2.4) 



The linearized normal fluid stress at the interface is, 

Pnn= lim /9[(9t¥3 + g?7(r,i)] ~ F^^, (2.5) 

where we have assumed the absence of externally imposed stream flow and that the only 
disturbances are in the form of dynamic surface waves. The fluid stress P„„ is balanced 
by spatially periodic varying material restoring forces from the surface scatterers, such 
as those due to surface bending and stretching. In the absence of fluid viscous stresses, 
the application of surface stresses must be carefully applied to ensure conservation of 
z— component forces, torques, and bending moments. These forces and the energy are 
manifestly conserved when the surface stresses are derived from an energy functional 
for the bulk fluid plus interface. For example Pzz and gradients of a{r) have to be 
considered in inviscid fluids to conserve wave energy. Thus, the net z-component surface 
stress balance reads 

[V^-(a(r-)Vi) - Vi(i?(f)Vi)] r;(f,t) = P., (2.6) 
which upon taking the time derivative and using (2.3) yields 



lim [Vx-(cT(f)Vi)-Vi(I?(f)Vi)]a,(^(r,z,t) = atF,,, (2.7) 

where D{r) and cr{r) are the flexural rigidity and surface tension respectively. In Eqn. 
( ^.7[ ), we have neglected the rotational inertial terms in the surface bending stresses. 
Recall that in thin plate theory, D = EdP — s^) where E and s are Young's modulus 
and Poisson's ratio of the surface material in its hulk phase. For open water between 



6 



Tom Chou 



the thin plates, Di — 0. A physical example is a large sheet of sea ice with alternating 
thickness. 

We will only consider, for all dynamical variables, time variations of the form e~*'^*. 
The neglect of rotational inertia terms is thus valid when uP' <^ Ep~^\~^, which holds 
for most applications. In the frequency domain, we combine the time derivative of the 
z— component stresses and use (2.3) to obtain 



lim [pu}^ip{f,t)-(pg-V^-{cr{r)Vi_) + V'i{D{r)V'i))dMr,z,t)]=0 (2., 



Equation (2^) and ( p.2| ) determine the velocity potential which has a e time depen- 
dence. The effects of spatially varying surface properties are implicit in the boundary con- 
dition (p.^). In the limit of uniform a{'r) = a and D{r, t) = D, ip (x e^''^ '" coshfc(/i -I- z), 
from which we obtain the standard dispersion relation 



= ( gk + -k^ + —kA tanhfc/i. (2.9) 
\ P P J 

This expression is valid far from (at least many wavelengths from) spatial inhomogeneities 
of the surface parameters (j{r) and 13 (r). 

However, when a and/or D are not uniform, surface waves diffract /refract from the 
regions of varying surface properties and the Fourier modes of the velocity potential at 
the interface mix with those of the surface variations. The rest of this paper deals with 
periodic variations in cr(r) and D{r) where the boundary condition ( ^.^ ) is to be used to 
solve V^<p. 

2.2. Block Functions and Periodic Solutions 

The following analysis is similar to the treatment of a single particle electronic wave 
function in a periodic potential. Details are given in Ashcroft & Mermin (1976). Another 
application can be found in photonic crystals where light propagation through a medium 
with periodic inclusions of different dielectric constants is studied, Joanopoulos, Meade, 
& Winn, (1995). These references describe in detail the terminology used; however, 
our problem contain some minor differences and we give a largely complete formulation 
below. 

Near sources or variations in a or D, the velocity potential must be decomposed into 
a complete set of eigenfunctions. The choice 



(p{r, z) 



^.p coshq{h -\- z) 



coshq/i 



(2.10) 



with q = lime^o Vy+^ manifestly satisfies Laplace's equation and (2.4). For the 
surface scatterers we consider one and two-dimensional periodicities, with lattice vectors 
a (see Fig. |l|). Since the surface properties are periodic in a, {i.e, D{r + a) — D{r)) they 
can be Fourier decomposed according to 



iG-r 



G 

D{r)^Y.DiG)e^' 



(2.11) 
(2.12) 
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where G are the corresponding reciprocal lattice vectors. The reciprocal lattice vectors 
for stripe, square, and triangular lattices with periodicity a are 



G 



27rTO 



27rm ^ 
G — X - 



47rm , 
G — ^ X 



mGo 



2?™ 



27rn 



(stripes) 



(square) 



(2.13) 



H — -^{x—Viy) (triangular) 



with m, n integers. Any particular reciprocal lattice vector in the set {G} for each 
type of lattice can be constructed by a countable number of others. We now exploit 
the periodicity of D{r) in the boundary condition ( ^.8|) and use Bloch's theorem which 
states that a function with periodic interactions and boundary conditions can be written 
in terms of a phase factor times a function invariant under discrete translations of the 
periodicity, if{f, z) = e^'^''^ 4>{r) . The function (l){r + na) = (/)(r) can be Fourier represented 
in a manner identical to D{r) or cr(r). Bloch's theorem and periodic boundary conditions 
on the entire array of scatterers (Born- von Karman boundary conditions) require q to be 
real and have the form (Ashcroft & Mermin (1976)) 



^Go(^) 



(2.14) 



where Ni is the large number of scatterers in the i*'* lattice direction, and Go{i) is the unit 
reciprocal lattice vector in the z*'' direction. The granularity of q is thus determined by 
the total number of scatterers that we apply the periodic boundary conditions to. Direct 
substitution of either e*« ''0(r) or ( |2.10| ) and ( ^.12D into (^7 
Fourier component 



yields for each independent 



where 



{^'-ni{G))^,-{G)= 

G'^G 

{q-G')-{q-G)\q-G'\a{G-G') 



q^G\^\q-G'\''DiG-G') + 
tanh \q-G'\hipi^{G') 



(2.15) 



Do, 



niiG) ^ g\q- G\ + -W- G\' + -^\q- Cr tanh \q- G\h 



(2.16) 



and (To = "■(0) and Dq 



P P 

D(0) are the spatially averaged surface properties. Equation 
(2.16) is the fundamental matrix equation that determines the relationship between fre- 
quency Lu and wavevector q. For the sake of mathematical clarity and simplicity, we 
will consider the case where a = relevant for long gravity waves (A > 2cm), or short 
bending waves (A < 2TT{g/Dy^'^). Thus the scattering originates only from the the peri- 
odically varying bending rigidities at the interface. The extension to include (t(G) ^ 
for studying surface waves in th e presence of periodic surfactant coverage is straightfor- 
ward. Note that in Eqn. ( 2.16 ), fq{G) is invariant to shifts in G. In other words there 
is only a set of Hi that uniquely determines ip^. Thus, the HiNi coefficients are 

coupled (</'g(0), 'PqiG), ip^{G'), etc.) by as many equations. We have therefore shifted 
the wavevector by G such that q — G falls within the boundaries delineated by the first 
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Bragg planes, i.e. the first Brillouin zone. Only the J^^ Ni number of (f s need be consid- 
ered, and we will eventually represent this by plotting in the reduced zone scheme, where 
q only takes on values in that part of the first Brillouin zone irreducible by symmetry. 
For example, in one dimension, the Bragg planes are points at — ±7r; however, we 
need only plot uo{qx) in the region < g^; < tt. 

Note that we are considering a linearly stable nondissipative system and all eigenval- 



ues poj"^ are real. The matrix (2. If) can be written in an inherently symmetric form by 
considering normalized Fourier coefficients f'^iG) = \q — G|tanh|g — G\h(p^{G). The 
problem is therefore a generalized eigenvalue problem A'lp'- — pw^Btp^ where B is diago- 
nal. Applying the inverse of the square root of B to both sides (Cholesky decomposition), 
the problem is recast into a symmetric eigenvalue problem, 

det\A - p{uj^ -nl)l\^0, (2.17) 

where 



A(G ^ G') = D{G - G')\q^ G\^/^\q- Gf tanh^/^ G|/i tanh^/^ G'\h. (2.18) 

The above symmetrization facilitates the calculation and computational convergence of 
the eigenvalues poj^. 

We now consider example one and two-dimensional periodic "surface potentials" D(G— 
G'). In one dimension, we consider square notch scatterers corresponding to a surface 
sheet with alternating bending rigidities D{x) = Di or D{x) = D2 depending on whether 
X is inside or outside a strip, respectively (see Fig 0(a)). For two-dimensional periodic 
surface scatterers, we consider both square (Fig|^(b)) and triangular (not shown) lattices 
of circular disks of radius i?o- Inside each disk, D = D2, while outside D — Di. The 
Fourier transforms in reciprocal lattice vector space for the above scattering potentials 
are, D{0) = £> = (!- f)Di + /D2, and 

D{G ^ 0) = sin(|G|a/) ID periodicity 

ttIGI 

(2.19) 

D{G^{)) = 2/ Ji ( I G I i?o ) 2D periodicity, 

\G\Ro 

where / is the fraction of area covered by the D2 material, and G is the appropriate 
reciprocal lattice vector for the surface periodicity under consideration. Our choices above 
are numerically the most laborious due to Gibb's phenomenon at the sharp discontinuities 



in D{f); nonetheless, we find the eigensolutions to (2.18) converge rapidly as a function 
of matrix size. 

As a consequence of a periodic differential operator, V\{D{r)V'\) , the band gaps in 
the water wave dispersion relation will, unlike periodic scattering in electronic, optical, 
and acoustic systems, more likely increase initially as frequency increases. This prediction 
is apparent when we consider relatively weak scatterers D such that ^IGji^ S> -D(|G|) 
where for certain G we can obtain an approximate analytic approximation for the band 
gaps in uP' . Using degenerate perturbation theory ai q — G near Bragg planes is tan- 
tamount to truncating A into block m x m matrices for m frequencies that are closely 
spaced. For simplicity, we consider one dimensional square stripes where doubly degen- 
erate frequencies near the n*'* Bragg planes {qx = 0, tt) are split according to 
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— jtanh(^ — hj (2.20) 

The stop band splittings increase with order as (I?2 — i^i)/!^ sin(n7r/), ahhough the 
absolute frequencies where they are found also increase as Drt' in this weakly scattered 
limit. Similar expressions can be found for the band gaps near Bragg planes in two- 
dimensionally periodic systems. 



D /n7r\4 ^ 



^ SI 



a / 



\ a / 



pnii 



sin(n7r/) 



3. Results and Discussion 



We numerically solve the eigenvalue problem represented by Eqns. ( 2.17 ) and ( 2.18 ) 
using standard methods (Press et ai, 1992). Typically, 256 plane waves (different wecG's) 
are taken for both the one and two-dimensional cases, such that the lowest 4-6 eigenvalues 
do not change appreciably (< 1%) upon doubling the matrix size. The calculations con- 
verge very efhciently for the lowest 4-6 eigenvalues especially in two dimensions because 
the extra Bessel function Ji make the off-diagonal elements fall off faster. The matrix 
calculations for smoother variations in D{f) become extremely simple; for example, A is 
tridiagonal when D is sinusoidal. 

Analysis of the following results will give a consistent qualitative picture of multiple 
scattering effects and surface wave propagation in the presence of periodic bendable 
plates. In our numerical plots, all distances are scaled with respect to lattice spacing, 
and wavevectors are measured in units of a~^. The frequency is nondimensionalized and 

4 



measured in units of ^ gja^ and the bending rigidity is measured in units of pga jix . 

3.1. One Dimensional Periodicity 

The choice Di = 0, D2 is appropriate for open water punctuated by flexible strips 
and illustrates the difference between gravity and bending waves. For concreteness, in the 
scattering strips, we choose D2 = 3.0, which for the Young's modulus and Poisson's ratio 
of sea ice {E = 6 x 10^°dynes/cm and s — 0.3), and periodicities of 1.0m, corresponds 
to a d ~ 0.8cm cover of ice. At this coverage, the dispersion within the ice sheet is 
predominately determined by bending forces (compare the terms gq and Diq^ / [pa^)) 
except at large wavelengths {q < n/2). 

The first few bands in the dispersion relation for wave propagating in an array of 
one-dimensional periodic surface scatterers are plotted in Figure ^j. Here, the depth 
h = 00, the filling fraction / = 0.6, and D2 — 3.0. The right panel shows the dispersion 
relation for gj, = as a function of Qx, the wavevector in the direction of periodicity. 
We have plotted the axis according to the reduced zone scheme were the dispersion 
relations are folded at every half unit reciprocal lattice vector Go, i.e. every Bragg 
plane. The band indices n correspond to the number of half reciprocal lattice base 
vectors Go/2 — xir/a included in the x-direction wavevector (the notation for q^ and 
the n*^ order band gaps is shown in Fig. ^). Note the wider higher order gaps in Fig. 
^. These gaps are simply discontinuities in the spectrum as wavevector parameter is 
increased and are associated with Bragg scattering planes. Physically, they are standing 
waves where the lower frequency branch corresponds to oscillations with most of the wave 
nodes(antinodes) located in the low(high) D{r) regions, Di and D2 here respectively. 
The higher frequencies are clearly associated with a translation of these low frequency 
standing waves, and have more nodes in the D2 regions. The left panel shows w as a 
function of Qy. However, in the reduced zone scheme, only the lowest branch (n = 0) 
corresponds to a wave propagating solely in the y direction. Interestingly, the motions 
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Dj=0.0, 0^=3.0, f=0.6 



CO 




Figure 2. The dispersion relation, or band structure of surface waves (with h — oo) in the 
presence of stripes of discontinuous D{x). Lengths are measured in units of periodicity a and 

frequencies are measured in units of \fgja. Surface values used are cr = 0, Di = 0, -D2 = 3.0, 
and filling fraction wja = f — 0.6. The right panel plots 0^(53;, 0), the left, (jj{0,qy), where 
a;(g) — uj{—q). The eigenfrequencies are plotted in the reduced zone scheme, where different 
bands n correspond to wavevectors displaced by nGo/2 — mi /ax. Note the first small band gap 
between tj„=o(7r,0) ~ 1.77 and aji(7r,0) ~ 1.79 



corresponding to this mode has no variation in the x direction. Note that this wave has a 
predominantly gravity wave-hke dispersion. The entire set of surfaces can be constructed 
by folding the page along the a; = axis and bending the left panel upwards. The surfaces 
are then bounded by the plotted curves. As a consequence of inversion symmetry, here, 
as well as in the two dimensional lattices, the dispersion relations are symmetric with 
respect to g — > —q. There are an infinite series of gaps in the qx direction, however, since 
the bands increase indefinitely in the qy direction, for any real uj, propagating modes in 
some direction can be found and complete {i.e, in all directions) band gaps do not exist. 
Complete gaps may exist in higher dimensions however, as they do in periodic optical 
and acoustic systems. 

Shown in more detail in Fig. |3| is the first band gap at qy = 0. The frequency range 
between ujoiT^TO) and cji(7r, 0) (labeled for Di = 0.03,1)2 = 3.0 only) correspond to 
frequencies where modes cannot propagate in the q^ direction. This fact can be shown by 
explicitly constructing the standing wave velocity potential within the gap frequencies, 
but can be seen more directly by considering the mode group velocities Cg = d^{q). 
Note that Cg^x ^ as a Bragg plane is approached due to the inherent symmetry of 
the problem. The sensitivity to Di 7^ is also apparent. Notice that for Di > 0, the 
dispersion relation develops an upward curvature indicative of increasing bending wave 
characteristics, especially in the bottom panel. This trend exists in a;o(0, qy) as well. 

We now explore the parameter dependence of the band structure. Figures ^(a) and 
(b) show the positions and widths of the first three band gaps as a function of depth 
h for fiUing fractions of / = 0.8 and / = 0.97 respectively. In both plots Di = and 
D2 = 3.0. The effects of h become important only in shallow water when \q — G\h < 1 and 




Figure 3. An expanded plot of the first band gap for various Di. (a). Solid curves, Di = 0.0; 
dashed curves, Di = 0.03. The repeated zone notation is indicated, (b). Solid curves, Di = 2.0; 
dashed curves, Di — 6.0 

manifests itself primarily through changing the excited wavelengths at fixed u. Therefore, 
the higher order gaps (higher \q — G\) saturate to their /i — > oo limit at smaller depths. 

The gaps are most sensitive to the filling fraction / and Di. Figure ||(a) shows the 
filling fraction dependence of the first three gaps for Di = 0.0 , D2 = 3.0, and h = 00. The 
lower plot. Figure ||(b), shows the same three bands when Di ~ 1.0. Band gap widths 
are not monotonic functions of / and generally increase near larger values of / when 
D2 > Di. As expected, the gaps are spread to lower / as Di is increased and actually 
are larger for smaller / when Di > D2 (not shown). All gaps disappear at / = 0.0 
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Figure 4. The first three band gaps as a function of water depth for 
Di = 0.0 , I>2 = 3.0 , ft = oo at (a). / = 0.8 and (b). / = 0.97. 



and / = 1.0 as the interface becomes uniform. Note that for the second (and higher) 
bands, the gap widths alternate and disappear for special values of /. This property is a 
consequence of the sinusoidal behavior of D{G). Higher bands have qualitatively similar 
behavior but are at much higher frequencies, beyond the scales of Fig ^. 

The variation in band gaps as a function of Di is also analyzed. The first, second, 
and third gaps as functions of D1/D2 {D2 = 3.0) for /i = 00 and / — 0.6 and / = 0.8 
are plotted in Figures ^(a), (b) and (c), respectively. The first gap, shown in Fig. ^(a) 
is most sensitive to variations in Di near Di = 0. The gap increases rapidly from a 
small constant value as Di is increased until it vanishes as Di approaches D2, when the 



surface becomes uniform. At this point the dispersion relation is governed by Eqn. 2.9 
in all directions. Figures ^(b), (c) however, show interesting reentrant behavior as Di 
is varied. There can be special values of Di 7^ D2 where the second and higher gaps 
vanish, similar to the vanishing of gaps as a function of / shown in Fig. |^. 
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Figure 5. The first three band gaps, u)i{-k,0) — (jJo(7I", 0) , tJi(0, 0) — tJo(0,0), and 
0^2(71,0) — iJi{n ,0), as a function of filling fraction / = w/a for (a). Di = 0.0 and (b). Di = 1.0, 
all other parameters are those of Fig H. 



3.2. Two Dimensional Periodicity 

The two dimensional dispersion relations depend in a more complicated way on the 
direction of q. We show the results for square and triangular lattices, where |Go| = 27r/a 

An/^/Sa and i?o 



and i?o — a^/Z/TT for square lattices, and |Go| = An/^/Sa and Ro = S^^'^a^J f /2'k for 
triangular lattices. The first lowest bands for waves propagating in square and triangular 
lattices are plotted along the paths in the Brillouin zone depicted in the corresponding 
insets in Figure (a) and (b) respectively. In the two-dimensional plots, we have chosen 
D2 = 100.0 and Di = 0.0, corresponding to d ^ 2.6cm of ice when a — IGOcm. We 
used / — 0.3 for the square lattice and / = 0.4 for the triangular lattice. Note that 
the maximum filling corresponding to close packing is / = 7r/4 and tt/2^/3 for the 
square and triangular lattices respectively. The dispersion relations shown in Fig. |^ may 
qualitatively represent monochromatic wave propagation through pancake ice zones. In 
Figure 0(a), the dashed curve corresponds to wave propagation over a uniform fiexural 
rigidity of D in the repeated zone scheme. The qualitative features depicted in Fig. |^ are 
not as sensitive to the parameters as the one-dimensional periodicity considered above, 
although oscillations in the band gaps at Bragg points X and M also occur, though are 
less pronounced compared to Lu{qx) in the ID model as / and D1/D2 are tuned. 

3.3. Complete Gaps, Boundaries, and Defects 

After a parameter search in /, Di , I?2, and h, we find no complete low order band gaps 
in two dimensions: There are always certain directions in which waves can propagate. 
This property is in contrast to photonic and phononic band structure, where complete 
gaps have been found by Joanopoulos, Meade fc Winn (1995)] , Kushwaha et al. (1993) 
and [Sigalas fc Economu (1992)| . The origins of this difference lie in nature of the wave 



scattering and the boundary conditions at the scatterer edges. Firstly, the dispersion 
caused by wave impedance mismatch due to variations in D is weak, i.e, since uP' ~ 
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Figure 6. The first three band gaps, (a), (b), and (c), respectively as a function bending 
rigidity mismatch D1/D2 at infinite depth and fiUing fractions / = 0.6 and 0.8. 



DiQ^ , Di must change by a considerable amount for the wavelengths 27: /q that would 
be excited in an infinite domain of uniform Di to change appreciably at any given uj. 
Secondly, although higher derivatives r}{r, t) are required to be discontinuous, the lower 
derivatives are continuous which limits the reflection from an edge discontinuity in D. 
Thus, scattering is not particularly strong in this system. Furthermore, the evanescent 
modes can couple among the scatterers for high filling fractions. This is a direct analogy 
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Figure 7. Two dimensional infinite depth band structure for (a), square and (b). triangular 
lattice of circular plates of radius _Ro- D2 — 100.0 , Di = 0. The filling fractions / used are 
indicated. 



to prism coupling of a totally internally reflected light ray at the edge of a prism. In 
the water scattering problem, these evanescent modes decay exponentially with a length 
scale set by the depth, (see Eqn. (A3)) except in the h ^ 00 limit where they decay as 
a power law. Thus, evanescent modes in infinitely deep water can especially be potent in 
effectively coupling the surface scatterings and prevent wave localization and complete 
band gaps. Since the gap widths in any direction tend to increase at higher bands, 
one might expect that complete gaps can be found at higher frequencies (or energies); 
however, this is not likely either since the dispersion relations and the positions of the 
gaps are very sensitive to — G at high frequencies, and the chances that gaps in all 
directions fall within each other are also diminished at high to. 
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An interesting correspondence exists for gravity waves (uniform interface) in the pres- 
ence of periodic bottom scatterers. The shahow water, depth averaged equation 



(3.1) 



is isomorphic to the acoustic wave equations in a fluid (no shear modulus). However, 
although complete acoustic band gaps have been found in square lattice geometries 
( ^igalas fc Economu (1992) ), Eqn. (3.1) actually corresponds to uniform compressive 
moduli, but different material densities in the fluid acoustic scattering context. In 
other words, the speed of sound and the density difference are constrained such that 
(ci/c2)^ = P2//O1 for the analogy with Eqn. (3.1) to hold. Under these conditions. 



3igalas & Economu (1992) suggest there are no complete band gaps. Therefore, we con- 
clude that the periodic bottom, shallow water wave problem has no complete band gap. 
A numerical search using the analogous matrix ( 2.16D for Eqn. (3.1) confirms the absence 
of a complete band gap. 

We have only quantitatively treated an infinite (periodic boundary condition) perfectly 
periodic array of scatterers. However, in real physical situations, boundary and defects 
in the periodicity exist. For example, waves generated in the open ocean may impinge 
on a field of periodic scatterers, or, defects may exist within a large array. Motions 
of scatterers induced by waves themselves can cause deviations from perfect periodic- 
ity, although this effect vanishes at the band gaps where only standing and evanescent 
modes are excited. Wave structure near edge and defect effects are formally described 
by imaginary wavevectors for frequencies within the band gaps, i.e, there are localized 
evanescent modes near boundaries and defects. In other words, the full dispersion relation 
is continued to the complex g— plane. These boundary and defects lead to localization 
phenomenon prevalent in many acoustic, electronic, and optical systems. 



4. Summary and Conclusions 

In this paper, we have presented a simple way of analyzing the dispersion relation 
for surface fiexural-gravity waves in the presence of flexible surface scatterers. The 



method is adapted from establishe d te chniques found in [Ashcroft fc Mermin (1976) 



Joanopoulos, Meade fc Winn (1995) , and Sigalas fc Economu (1992) for calculating mul- 
tiple scattering in electronic, optical, and acoustic phenomena. Figures encapsulate 
the main conclusions of this paper. Our calculations have implicitly included all scat- 
tered evanescent modes by solving the problem in discrete Fourier space, and we have 
treated the most computationally difficult case, that of sharp surface discontinuities. We 
conclude that: (1). Band gaps can exist in flexural surface scattering of monochromatics 
waves, (2). The band gaps generally tend to increase with excitation frequency uj, (3). 
The sensitivity do depth is weak, (4). The sensitivity to / and D1/D2 can be strong, 
where (5). Certain special values of / and D1/D2 have vanishing gaps, and (6). The 
high dispersion in this system prevents the formation of complete band gaps in the two- 
dimensional case. The gaps discussed discussed here are associated with nonpropagating 
waves and correspond to the "fully resonant" scattering from a flnite number of periodic 
bottom undulations discussed in Davies, Guazzelli & Belzons 1989 where their analyses 
broke down. Our approach is inherently applicable in the fully reflecting, or resonant 
regimes (band gaps) although for an infinitely periodic system of flexible surface scatter- 
ers. Further exploiting the acoustic analogy, we also conclude that complete band gaps 
do not appear in periodic bottom, shallow water wave scattering. In Appendix A, we 
related our results with the calculations of reflection R, and transmission T coefficients 
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of scattering from a finite number of surface flexible obstacles. In particular, the rela- 
tionship between |r(ti')| and its phase 6{uj) was found at all frequencies outside the band 
gaps. In Appendix B, we considered the possibility of a uniform flow underlying the 
periodically modulated ice sheet. 

Our description of band gaps, frequencies where water waves cease to propagate, imply 
possible structures for wave damping applications such as breakwaters. Within a large 
field of periodically spaced flexible plates, the band gaps determine the standing wave 
response, and reflects all travelling waves from the large fleld of plates. If a structure 
embedded in this fleld is susceptible to transverse wave forces, then it may be desirable 
to have large band gaps where modes are standing. However, if the periodic surface 
field is to function as a wave frequency filter, then band gaps are to be avoided in the 
region of structural resonant frequencies. The important result (5). above suggests that 
band gaps are generally not monotonic in / and D2/D1 and there are special values 
of these parameters which are to be targeted or avoided depending on the application. 
Consideration of whether wave generation (by e.g, wind) occurs inside the array should 
be taken as well. Although it appears that no simple complete gaps exist for periodically 
bottom or surface scattered waves (due to a large dispersion in the surface scattering 
case, and due to uniform gravity in the shallow bottom scattering case), by appropriately 
choosing the two-dimensional surface periodicity relative to mean wind directions, one 
may nonetheless shunt propagating waves for desired band gap frequencies. The effects 
of directionality of externally impinging and internally generated waves can be inferred 
from the dispersion relations plotted in Figures || and 0. More complicated structures 
where the basic unit cell of periodicity contains more that one scatterer can be treated 



by similar methods as described in this paper and elsewhere (Ashcroft & Mermin (1976) 
and Joanopoulos, Meade fc Winn (1995)] ) 



Applications to wave propagation in Marginal Ice Zones (MIZ) may also be in- 
terpreted in terms of multiple scattering (see for example [Fox fc Squire (1990) and 



3quire et al. (1995)), with the perfect periodicity treated in this paper as an extreme 
limiting case. Defects in the perfect periodicity can give rise to evanescent local- 
ized modes which may be relevant for "Roll-over" and other anomalous effects seen 



in measurements of wave energy in the MIZ as discussed by Squire et al. (1995) and 



Wadhams et al. (1988). Further work may extend the analyses to specific problems in- 
volving wave forces on surface ice, ice breaking, and wave channelling. It may be feasible 
to break ice in certain patterns to channel surface water wave energy much like optical 



waveguides ( Joanopoulos, Meade fc Winn (1995) ) 



The author thanks TJ Pedley for bringing to his attention the work of DV Evans, MC 
Payne for helpful discussions, MG Worster for comments on the manuscript, and Dirk 
Laurie for checking the numerical code. The author also acknowledges the Wellcome 
Trust for their support. 



Appendix A. Reflection and Transmission in ID 

In the one dimensional periodic scattering problem, the existence of band gaps can be 
explicitly calculated from considering the reflection R, and transmission T coefficients 
from a finite number n of scatterers. The calculation of R and T in general water wave 
scattering systems is mathematically difficult, however, in all cases, two of the four un- 
knowns are eliminated via conservation of energy and wave momentum. Here, we derive 
a relationship between the scattering coefficients and phases of scattering from a finite 
number of scatterers to the exact periodic solutions computed above, and show that in 
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Figure 8. N orm al incidence scattering from a single (n = 1) surface scatterer. The 
relationship Eqn. (A 9) is valid for any general symmetric one-dimensional scattering problem. 



certain reasonable limits, a third condition arises which allows efficient calculation of the 
reflection and transmission amplitudes and their phases, given any one of them. There- 
fore, this relationship may be useful when combined with variational calculations of the 
reflection and transmission amplitudes. This third constraint may also be useful exper- 
imentally, when only one of the four unknowns (for single one-dimensional scattering) 
can be readily measured. 

A single scatterer of general shape D{x) is depicted in Figure |^ with the only provisions 
that D{\x\ > w/2) — Di is constant and D{x) = D{—x). Consider n replicas of such 
surface scatterers centered about a; = 0. Note that the most general solution to Laplace's 
equation for |a;| > (n — l)a/2 -|- w/2 can be written as a complete set of eigenfunctions 



(f{x, z, bj) = A±e^^'^"'' cosh fco(/i + + ^ B^e^'^"""' cos Hn{h + z) (A 1) 
where fco is the positive real root of 



n=0 



/cq -I- (7 fco I tanh k^h 



and the defined by 



Dl 5 



gKn I tan(K„/i) — uj 



(A 2) 



(A3) 



Thus, the far field (|a;| oo) asymptotic solution consists of travelling (again, the e~"^* 
factor has been suppressed) propagating waves, with the coefficients A± related to R and 
T. In a single scattering problem, we are not interested in the complicated details of the 
evanescent components -B^ , except for the fact that they modify the potential and match 
near the scattering region ^ {n— l)a/2 + w /2. Since nonpropagating modes decay 
exponentially according to Eq. ( |A 3| ) , the longest-ranged mode decays exponentially with 
a length scale comparable to the depth h. For infinite depth fluids, the infinite number 
of poles in the complex k plane collapse into a branch cut with the result that the 
nonpropagating mode now decays ~ (fco|a;|)~^. We assume that the wavelength 27r/fco or 
h is small compared to a — w which can be achieved by waves excited at high frequencies 
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and/or considering small /. The velocity potential in a periodically repeated structure, 
in the region \x\ < na/2 can be thus written as combinations of single scattering problems 
with incident waves from x < —na/2 and x > +na/2, where 



- — p^Ko 



and 



Te 



fcox coshfco(fe + z) ^ 
cosh kgh 

,koxcoshko{h + z) 



cosh fco(/i + z) ^ , . ^ , , ,„ 

j/cox — X < -(n - l/2)a - w/2 



cosh feo^ 



cosh koh 



x>{n- l/2)a + w/2 
(A 4) 



^+ ^ ^-»fco:. coshfco(/i + z) 



Te" 



cosh fco/i 
fe^^cosh fco(/i + 2:) 



ikox 



coshfco(/i + 
cosh fco ft. 



+ Ci(-a;,z) a; > (n- l/2)a + u;/2 



cosh koh 



+ ^2{x,z) x<-{n-l/2)a-w/2 

(A5) 

where the functions and ^2 denote the evanescent contributions which in general 
depend on R and T, decay away from the scatterer, and are small at a; = ±na/2. The 
entire potential in the periodic system is thus 



ip{x,z) ^ Ci(p- +C2ip^. (A 6) 

We apply Bloch's theorem to the potential and Vx for all z across n scatterers and obtain 



(p{-na/2) = e'"«-"^(na/2) 

(A 7) 

dxV{-na/2) = e'"«-'^a^^(na/2). 

The ^1^2 are chosen such that ( |A 7| ) determine Ci and C2 for all z. Considering a region 
straddling the n scatterers far away from the nonpropagating modes; conservation of 
energy and momentum across its surface (a; ~ ±na/2) requires 



|i?|2 + |T|2 ^ 1, 

1- R = T 



{Al 



respectively. Defining the complex values as i? = |i?|e'-*^ and T = |T|e**, (A 8) immedi- 
ately give X = ^± 7r/2. Upon substitution of (A 6) into (A 7) and using (A 8), x, and S 
to reduce the result, 



cos(nfcoa + S) 



m 



= cosnq^a + 0(^i,2(na/2)) + 0(^1.2)- 



(A 9) 



Equation ( |A 9] ) represents a con ditio n between |r| and 5 when fco(w) and Qx are given. 
One can readily obtain fco from (A 2) and from the results in Section 3. However, only 
higher bands (higher frequencies) should be used, as these are associated with sufficiently 
high wavevectors to allow neglect of the evanescent ^1^2 terms away from the scattering 
edge. Consider the reduced frequency u) ~ 5.45, and the reduced qx = n (in the repeated 
zone scheme) of the fifth band (therefore wavelengths A < a/5) in Fig. 0. We have 
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cos (nfco(5.45) + (5(5.45)) ~ (-l)"|r(5.45)|, (A 10) 

where ko is to be found from (A 2). In general, for any fixed qx, there are an infinite 
set (neglecting the lower bands) eigenfrequencies lo such the the functional relationship 
between |r(ix')| and (S(w) are identical. The relation (A 9) holds for longer wavelengths if 
h <^ a since the evanescent modes are shorter ranged. 



Appendix B. Uniform Flow 

The analyses given in this paper can be easily generalized to cases with uniform current 
underneath the surface scatterers. The related bottom scattering problem has been 
studied by Kirby, 1986. The linearized interfacial normal stress due to the fluid in this 
case is 



Pr,n = p\dt^p + Uf)-V^Lp + gr]y (Bl) 

where Uq is the uniform stream velocity just beneath the surface. Further derivation will 
be simplified by considering only a fully material-covered interface, where the surface 
deformations follow Lagrangian coordinate and are not convected by Uq. The kinematic 
condition in this case is identical to Eqn. (|]J). Thus, for the sake of simplicity, we do 
not consider open water, but can allow for varying thickness in ice cover for example. 
Since we are considering an ideal fluid, the shear layer becomes infinitely thin and Uq 
contributes only to the normal stress. The effects when there is open water can be easily 
included by considering Uo{G) as nonzero constant in the open water regions, and zero 
underneath the scatterers in the kinematic boundary condition only. For a fully plate- 
covered interface {a — 0), the additional term Uq in the z component stress balance leads 
to the generalized eigenvalue equation 



det|A- [lo^ -ni{G)-ujUo-{q-G)j 1| = 0, (B2) 
which upon defining unp^{G) = ipq{G) can be rewritten in terms of a standard eigenvalue 



problem for {ipq'{G), ijj^{G)) as described in Press (1992). Clearly, only the component of 



q — G perpendicular to Uq will have inversion symmetry. The band structure will now 
be more complicated, and the band gaps will in general be different depending on the 
wavevector direction relative to Uq- 
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